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Abstract 

The Bray-Humayun model for phase ordering dynamics is solved 
numerically in one and two space dimensions with conserved and 
non conserved order parameter. The scaling properties are anal- 
ysed in detail finding the crossover from multiscaling to standard 
scaling in the conserved case. Both in the nonconserved case 
and in the conserved case when standard scaling holds the novel 
feature of an exponential tail in the scaling function is found. 
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1 - Introduction 



In this paper we are concerned with scahng behaviour in the phase or- 
dering dynamics of a system quenched below the critical point [1]. Specif- 
ically, we consider a system with an A'"-component order parameter = 

{(f)i{x), 0Ar(x)) quenched from high temperature to zero temperature whose 
dynamics are described by the zero noise Langevin equation 
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where p = for non conserved order parameter (NCOP), p = 2 for conserved 
order parameter (COP) and V{(f) = + is the local potential with 

{r < 0, g > 0). One of the reasons for the continuing interest in this type of 
problem is that a theoretical derivation of scaling on a first principles basis 
is still lacking except for a few exactly soluble models[2,3]. 

Let us first give a qualitative description of what goes on during the phase 
ordering process. Initially the system is prepared in a high temperature state 
where the order parameter is spatially uncorrelated 

< , O)0^(f' , 0) >= A5^fs(x - f ) (1.2) 

with a,P — 1,...,N and A is a constant. The local order parameter prob- 
ability distribution has as a peak of width A centered about (pa — 0. As 
the quench develops there is first a fast process (early stage) where this 
probability distribution, after a short time to, relaxes to equilibrium in the 
local potential, depleting the origin and developing a peak structure all 
around the bottom of the potential. In this span of time the local vari- 
ance < 4>l^{x,to) >= S{to) loses memory of the initial condition A reaching 
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a value very close to the final saturation value Seq = —r/g. At this point 
the system is almost at equilibrium, namely ordered, on a short length scale. 
The subsequent time evolution amounts to coarsening of the ordered regions 
in order to reduce the excess interfacial free energy. During this process 
(late stage) the only important time dependence is in the linear size of the 
ordered regions which typically grows according to a power law L{t) ~ t^^^ 
with 2; = 2 for NCOP and z = 3 or 2; = 4 for COP respectively with N — 1 or 
N > 1. When L{t) is large enough to dominate all other lengths the quan- 
tities of interest exhibit scaling. The main observables are the equal time 
order parameter correlation function G{f,t) =< (f)a{x,t)(f)a{x + f,t) > and 
the structure factor C{k,t) obtained by Fourier transforming G{f,t) with 
respect to space. The local variance of the order parameter is related to 
these quantities by S{t) = G(f = 0,i) = / j0^C{k,t). Actually, we will be 
interested in the quantity 

R{t)^r + gS{t)^g[S{t)-S,,] (1.3) 

which monitors how the saturation value of S{t) is reached. 

According to the scaling hypothesis, all time dependence can be expressed 
through L{t). The dominant behaviours for large L{t) are given by 

R{t) = ~ (1.4) 

with ^ = 2 for systems with vector order parameter [4] and 

C{k,t)^SeJ{L{t),kL{t)) (1.5) 

where the function f{L, kL) must go over to 5{k) in the limit L — > 00 in 
order to reproduce the Bragg peak corresponding to the final ordered state. 
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In other words f{L, kL) is a smoothed out 5-function on the scale L{t) with 
the normahzation property 

J -^j{m,kL{t)) = i. (1.6) 

Experiments, numerical simulations and soluble models, with the exception 
of the exact solution of the large- model with C0P[3], yield the following 
form of scaling 

f{L,kL) ^ L'^F{kL) (1.7) 

which we refer to as standard scaling. By contrast, when the model with 
N = oo and COP is solved one finds out that this pattern of scahng is not 

obeyed. In that case the system behaves differently since next to L{t) there is 
another divergent length k^}{t) ~ L/(lnL)^/^ where km{t) is the peak wave 
vector of the structure factor. Consequently, there appears a logarithmic 
correction in (1.4) 

i?(t) = --A_(inL)V2 (1.8) 
and in place of (1.7) one has the qualitatively different form 

fmk/k^)^{L'k'-'')^^'''-^ (1.9) 

with '4>{x) — l — {x^ — Vf. This pattern of scahng is referred to as multiscahng. 

It should be stressed that the above solution of the N — oo model is the 
only available analytical solution of a system with COP. Hence, due to the 
difficulty of discriminating between multiscahng and standard scaling on the 
basis of the usual data collapse analysis, one may reasonably ask the question 
whether multiscahng might in fact be a generic feature of all systems with 
COP. In other words, putting x = kL and neglecting logarithmic differences 
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between L and k^, i.e. letting L = u/km with u constant, one can postulate 
the general scaling form 

fmx) = [L{t)Y^-^F{x) (1.10) 

which contains (1.7) and (1.9) as particular cases respectively with (/?(x) = d 
and (p{x) = dip{x/u). It is then matter of computation or experiment to 
extract the spectrum of exponents ip{x) and to check whether it is fiat like in 
standard scaling or there is a genuine x-dependence implying multiscaling. 

This kind of analysis has been carried out on the data for the structure 
factor obtained from the simulation[5,6,7] of systems with COP and with N 
ranging from 1 to 4 in two and three dimensions. In all of these cases the 
observed behaviour of (fi{x) is consistent with the flat spectrum characteristic 
of standard scaling. Furthermore, analytical work of Bray and Humayun[8] 
(BH) on a model with N large but finite and d > 1 suggests that standard 
scahng holds for any finite N, while multiscaling is only a feature of the 
special case N = oo. Actually, the picture that BH put forward is that there 
exists a crossover time t* which depends on N and the initial condition A 
and such that multiscaling holds for t < t* whilst for t > t* standard scaling 
sets in. Therefore, different asymptotic behaviours are obtained according to 
the order of the limits t ^ oo and N ^ oo, as it was conjectured very early 
on by Yoshi Oono[9]. If the limit t — > oo is taken first, standard scaling is 
observed asymptotically for any value oi N, as long as N < oo. Conversely, 
if the limit — > oo is taken first, then the crossover time t* diverges and 
the asymptotic behaviour exhibits multiscaling since the regime of standard 
scaling can never be reached. 

What is at stake in this question of standard scaling vs. multiscaling is the 
nature of the symmetry underlying scaling behaviour[5]. The results of the 
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simulations could be regarded as non conclusive, since one could well imagine 
a spectrum (p{x) which is dependent on N and which interpolates between 
the N — oo behaviour and the standard scaling behaviour as N gets smaller 
and smaller. Then, for values of N of order unity, such as in the simulations, 
it might be difficult to decide whether a (p{x) with a weak dependence on x is 
evidence for standard scaling or for multiscaling. Instead, the result of BH is 
clearcut and states that the symmetry underlying asymptotic dynamics leads 
to standard scaling for any finite N. This result is quite important from the 
point of view of theory since theoretical progress in this field so far has heavily 
relied on the use of very clever but uncontrolled approximations [10]. This is 
due to the difficulty of developing systematic and controlled approximation 
schemes. An exception is the l/A'"-expansion for systems with NC0P[11]. 
The result of BH ehminates the possibihty of extending the l/A'"-expansion 
as a systematic expansion scheme to the conserved case. 

For the relevance of this issue, in this paper we have made a detailed study 
of the crossover from multiscaling to standard scaling through a comparative 
analysis of the numerical solution of the BH model with NCOP and with 
COP. Our aim is to proceed to an unbiased analysis of the scaling properties 
in order to have a check on the BH picture without any a priori assumption 
on the type of scaling, and to analyse in detail the difference between the 
conserved and non conserved case. In this respect our work is quite different 
from that of Rojas and Bray[12]. These authors do perfom a numerical 
solution of the BH model after the standard scaling ansatz has been made, 
while we first solve for the structure factor and then we proceed to the scaling 
analysis on the basis of the uncommitted general form (1.10). 

The paper is organized as follows: in section 2 we present the model and 
we elaborate on the difference between standard scaling and multiscaling. 
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introducing the observables best suited to distinguish one from the other. In 
section 3 we illustrate the method of solution with a test of its validity made 
by comparing numerical data with the analytical solution in the N — oo case. 
In section 4 we present the results for finite N in one and two dimensions 
and in section 5 we make some concluding remarks. 



2 - BH Model, Standard Scaling and Multi- 
scaling 



By using the gaussian auxiliary field method of Mazenko[13], BH have de- 
rived [8] from (1.1) a closed equation of motion for the equal time correlation 
function within the framework of the 1/iV-expansion. Retaining non linear 
terms up to first order in 1/iV one has 



dG{f, t) 
di 



2(iV)^ 



W'G - R{t) (G + -G 



(2.1) 



where R{t) is a function of time which must be determined by the equilibrium 
requirement 

lim ^(f = 0,t) = Seg = -r/g. (2.2) 

The corresponding equation of motion for the structure factor is obtained 
after Fourier transforming with respect to space variables 



dC{k,t) 
di 



= -2F [k^ + R{t)\ G{k, t) - 2—R{t)D{k, t) (2.3) 



where D{k,t) is the Fourier transform of G^{f.t). Notice that in the limit 
N —>■ oo (2.3) reduces to the equation which has been studied in [3] 

dG{k, t) 



dt 



= -2kP 



k^ + R{t) G{k,t). 



(2.4) 



In this latter case R{t) is defined self-consistently by (1.3). We shall retain 
this definition of R{t) also in the finite N case since from (2.3) follows that in 
order to reach equilibrium R{t) must vanish and with the definition (1.3) the 
condition \imt^^R{t) = is an implementation of the requirement (2.2). 

Although (2.1) or (2.3) have been derived by a truncation procedure based 
on the l/A'"-expansion, the solution of the equation is not of first order in 
1/A^, since it contains all orders in 1/N. Actually, it is not possible to assess 
precisely what is the relationship of this solution with the systematic 1/A^- 
expansion performed on the basic equation of motion (1.1). Presumably, it is 
some kind of infinite partial resummation intertwined with the uncontrolled 
approximation inherent in the use of the gaussian auxiliary field method of 
Mazenko[10,13]. Hence, (2.1) or (2.3) should be regarded as the definition of 
a model, the BH model, for phase ordering dynamics with an A'^-component 
vectorial order parameter which in the A" — > oo limit reproduces the usual 
large- A^ limit of (1.1) for the dynamics of the structure factor. 

In order to make the scahng analysis of the BH model, let us integrate 
formally (2.3) from some instant of time to onward 



where Q{t) — fQdt'R{t'). Choosing to in the scaling region, according to the 
standard scaling hypothesis we have the asymptotic behaviours 



C{k,t) = C'(^,io)e-'['='^'(*-*°)+*^'(^(*)-^(*°))] 




(2.5) 



C{k, t) 

R{t) 



SeqL\t)F{x) 
-hL-\t) 



(2.6) 

(2.7) 




, for NCOP 
, for COP 



(2.8) 
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with X = kL{t), L{t) = = L(to) and 6 is a positive constant to be 

determined. 

Let us first consider the case of NCOP. Performing the above ansatz on 
(2.5) we find 

F{x) = F(xo)(a;/a;o)^'-'^e-2("'-"o) 
r ^i^/^r-'e-''''-'''^V{x') (2.9) 

iV Jxo X 

where Xq = kLo and V{x') = J ^^F{\ x' - Xi \)F{\ Xi - x^ \)F{x2). 
The condition (2.2) requires 

j^.n^) = ^ (2.10) 

which gives an equation for h. Letting — > and requiring (2.10) to be 
satisfied, in the N — oo case we obtain Ab — d — and 

F{x) = F(0)e-2^' (2.11) 

while for finite N we find 46 — d < 0. 

Let us now go to the COP case. Making the scaling ansatz into (2.5) 
with p = 2 we find 

F(x) = F(xo)(xo/x)'e-^^^^"-^o)-M-'-4)] 

+ ^ r cia;V'^+^e-^[(-'-^")-^^(^^-^'^)]p(a;'). (2.12) 
Nx'^ Jxo V / V / 

Now, if we let again ^ 0, in the N — oo case F{x) vanishes identically and 
it is not possible anymore to satisfy (2.10). This is the breakdown of standard 
scahng in the large- limit which leads to multiscahng[3] . Conversely, if N 
is kept finite, (2.12) yields 

F{x) = ^ r dxV'^+^e-^[(^'-^")-^^(^^-^")]p(x'). (2.13) 
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This is the equation that BH have solved finding a non trivial solution and 
reaching the conclusion that for any finite N standard scaling holds also for 
systems with COP. 

According to the above discussion the first term in the right hand side of 
(2.5), the one which survives after — > oo, is responsible for multiscaling 
behaviour while the second one is responsible for standard scaling. The 
competition of these two terms is expected to generate a crossover time t* 
such that multiscahng behaviour of the type found with N — oo holds for 
t < t* while standard scaling eventually sets in for t > t*. 

In the following we will make a numerical study of the scaling properties 
of the structure factor in the BH model on the basis of the general scaling 
form (1.10). The primary interest is in the discrimination between standard 
scaling and multiscaling and in the study of the crossover. The analysis will 
be carried out through the behaviour of the spectrum of exponents (fi{x) as 
described in the Introduction and through the behaviour of the quantity 

Y{t) = -R{t)L^{t) (2.14) 

which discriminates between standard scaling and multiscaling on the basis 
of the asymptotic behaviours 

YU\ [ ~ H^) 1 standard scahng ,^ . 

^ ^ \ ^ (Int)^/^ , for multiscaling. ^ ' ' 
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3 - Method of solution and numerical results 
for N = oo 



In order to investigate the scaling properties of the BH model the discretized 
version of (2.1) is integrated numerically via a simple finite difference first 
order Euler scheme. The initial condition is given by G{f, 0) = A for r = 
and G{r, 0) = elsewhere. The boundary conditions are chosen to be 
periodic, but open conditions have been tested not to affect the final results. 
Prom the values of G{r,t) the structure factor is then obtained via Fast 
Fourier Transform and from these two functions all quantities of interest are 
computed. 

Two opposite requirements enter in the choice of the parameters of the 
numerical solution, and in particular of the linear dimension L of the system. 
A large number of lattice sites is desirable to avoid that the discretization 
of space may hide the subtle difference between standard scaling and mul- 
tiscaling. On the other hand, fewer sites speed up the computation and 
investigation of later times is feasible. Resorting to parallel computing we 
have managed to perform the numerical integration on large systems and for 
sufficiently long times. In principle, one can solve (2.3) to obtain directly the 
structure factor, but the presence in (2.3) of a double convolution integral, 
which cannot be parallelized efficiently, makes this alternative computational 
scheme much slower. 

For all runs the value of the mesh size has been taken Ax = 1, while the 
time step At has been changed depending on the values of and d in order 
to prevent numerical instabilities. In particular, for NCOP, At = 0.01 for all 
values of d and N except when N — 10. In this latter case, we have taken 
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At = 0.005. For COP, At = 0.05 for d = 1 and At = 0.01 for d = 2. For 
the parameters of the potential V{4>) we have chosen the values r = — 10 and 

— * 

After computing the structure factor C{k,t) for several different times, 
the spectrum of scaling exponents (p{x) can be obtained by using the general 
scaling form (1.10), which can be rewritten as 

In C{k, t) = cp{x) In L{t) + In F{x) (3.1) 

where L{t) = and x = kL{t). Hence, plotting \nC{x/ L{t),t) vs. InL(t) 
at fixed x one can measure ^p{x) from the slope and F{x) from the intercept 
with the vertical axis. However, from the numerical point of view it is more 
convenient to use a slightly different procedure, because x/L{t) could turn 
out to be too small or too big with respect to the available values of k. For the 
NCOP case, C{k,t) is computed not as C{x/ L{t),t) but as C{xk2{t),t) where 
k2{t) is defined by C{k2{t),t) = C{0,t)/2. This introduces in Eq. (3.1) an 
additional constant term given by the logarithm of the proportionality factor 
between L{t) and k2{t). Furthermore, the error in the determination of k2{t) 
and C{xk2{t),t) is greatly reduced by the use of linear interpolation between 
the discrete values of k. In the COP case C{k, t) is computed as C{xkm{t),t) 
where A'",„(t) is the peak wave vector of the structure factor. The logarithm of 
C{xkm{t),t) is plotted versus log{L'^{t)k'^'^{t)). With this choice, for N — oo, 
the slope (fi{x) is given by d%l){x) rather than by dip{x/u). This makes the 
comparison between numerical and analytical results easier also for A?" < oo. 
Again the quality of the fit is enhanced by determining the peak wave vector 
and C{xk„i{t),t) via cubic and linear interpolation, respectively. 

In order to check the quality of the numerical method, let us consider the 
N = oo case where exact analytical results are available. We solve for C{k, t) 
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in one and two space dimensions with A ranging in the interval (0.01, 10). 
We discuss first the case of NCOP and then the case with COP. The moti- 
vation for doing this computation is also to establish clearly the behaviour 
of observables according to standard scahng (NCOP) and to multiscaling 
(COP). 

NCOP 

In order to analyse the behaviour of Y{t) in Fig.l \nY{t) has been plotted 
versus In(lni) for d = 1 and d — 2. In both cases laY{t) displays the ap- 
proach to the asymptotic constant value lnd/4 of standard scaling predicted 
by (2.15) through a transient dependent on the initial condition A. No de- 
tectable dependence on the initial condition is found in the behaviour of ip{x) 
which in agreement with (2.6) follows the constant behavior ^p{x) = d. Simi- 
larly, the numerical results for the scaling function reproduce accurately the 
gaussian behavior (2.11). 

COP 

With COP the behaviour of \n.Y{t) is qualitatively different from what we 
had above with NCOP. In place of the relaxation to a constant value now 
(Fig. 2) there is an upward increasing trend revealing multiscaling. In the 
time of the computation there is still a dependence on the initial condition 
A, with a faster convergence to the asymptotic behaviour ~ ln(l/21n(t)) 
given by (2.15) for higher values of A. The asymptotic behaviour has not 
been reached in the time of the computation due to the much slower dynamics 
of COP. 

Multiscaling is most clearly illustrated by the behaviour of ^{x). It is 
interesting to see how the spectrum of exponents depends on the time interval 
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of observation. In Fig. 3 the evolution of (p{x) in subsequent time intervals 
has been plotted for different values of A and for d = 1. Results for d — 2 
are similar. Fig.3 demonstrates the relaxation of (f{x) to the asymptotic 
behaviour given by (f{x) — d%l){x). As remarked above the relaxation is 
faster for higher values of A. The late stage results are displayed in Fig.4 
both ior d— 1 and d — 2 showing the independence from A of the computed 
i^{x). This suggests that, at least in the range of x considered, ip{x) reaches 
the asymptotic regime faster than Y{t). 



4 - Results for finite 



In this section we illustrate the solution of the BH equation with finite N 
obtained by the numerical method described in the previous section. 

NCOP 

The standard scaling behaviour of systems with NCOP is manifested (Fig. 5) 
first of all in the behaviour of Y{t) which according to (2.15) goes to a con- 
stant asymptotic value h{N) smaller then d/A and decreasing monotonically 
with N . The transient preceding the asymptotic behaviour now depends 
both on A and N. Asymptotic behaviour independent of A and N instead 
is manifested by ip{x) which displays with great accuracy the flat behaviour 
i^{x) = d. A significant N dependence shows up (Fig.6) however in the scal- 
ing function F{x). According to the analysis of section 2, a deviation from 
gaussian behaviour is expected for finite N in the tail of F{x) due to the 
second term in the right hand side of (2.5) and this deviation clearly should 
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be more important for small values of N. The plot of In F{x) vs. x reveals 
the interesting feature that the tail decays exponentially rather than follow- 
ing the generalised Porod's law ~ x~^'^"'"^^[14]. Simulations of a system with 
NCOP and N > d have been performed by Toyoki[15] finding a tail which 
decays with a power much higher than that of the generalised Porod's law. 
Our result suggests that an exponential fit might be appropriate also in this 
case. 

COP 

The picture is more complex and interesting with COP. Fig.7 and Fig. 8 
display respectively the behaviour of lnY(t) for a fixed value of A'^ with 
varying A and viceversa for a fixed value of A with varying A^. What emerges 
from a comparison with the analogous data for = oo is that for fixed 
(here A^ = 300) the behaviour is of the standard scaling type for A sufficiently 
small (e.g A = 10~^) while it is of the multiscaling type for A large (A = 10) 
with an interpolating behaviour for intermediate values of A. Similarly, 
for A fixed (A = 0.01) the behaviour goes from standard scaling for A^ = 
300 toward multiscaling as N grows very large. This pattern fits with the 
crossover picture illustrated in the Introduction allowing for a crossover time 
t* which grows both with A^ and with A. Standard scaling then applies when 
the values of A^ and A arc such that t* is very short. For higher values of 
A^ and A, instead, t* can be made long enough for the system to develop 
multiscaling behaviour before the asymptotic standard scaling behaviour is 
reached. If only multiscaling behaviour is observed, as for instance for A'^ — 
300 and A = 10 or for N = 10^ and A = 0.01, it means that for those 
values of N and A the crossover time t* is larger than the maximum time 
reached in the numerical computation. From these data it is very difficult. 
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though, to infer the quantitative dependence of t* on N and A. BH have 
proposed[8] the analytical form t* ~ (AA^)'^/'^(lniV)'^, which holds for d > 1. 
However, the predictions from this formula seem to be quite off from what 
we observe. For instance for N — 300, A = 1 and d — 2 the above formula 
gives t* ~ 10^ which is large enough to expect an observable crossover, while 
for these values of the parameters we find only standard scaling behaviour 
(Fig.7 and Fig.9). 

As we have seen with N = oo the distinction between standard scaling and 
multiscaling is most effectively manifested through ^p{x). Thus, according to 
the crossover picture obtained from Y{t), it should be possible to produce 
standard scaling or multiscaling in ^p{x) by properly choosing the values of N 
and A. In Fig.9 we have analysed the evolution of (f{x) in time for A = 0.01 
and different values oi N ior d — 2 (similar results are obtained for d — 1). 
For N ranging from 10^ to 10^, (p{x) displays standard scahng over all time 
intervals implying that t* is of order one. For larger values N — 10^, 10^ one 
can definitely recognize multiscaling type of behaviour over the initial time 
intervals evolving toward standard scaling in the later time intervals. Here it 
is difficult to assess the value of t*, but it must be of the order of magnitude 
of the time of observation. Finally, for = 10'' the behaviour of ^p{x) is of 
the multiscaling type over all time intervals implying that t* is larger of the 
time of computation. 

In order to complete the analysis in Fig. 10 we have plotted the logarithm 
of the scahng function F{x) vs. x for different values of N finding again an 
exponential tail hke in the NCOP case. In this case there are small secondary 
peaks superimposed on the tail which scale like L{t) and which become more 
pronounced as N grows. Exponential tails have been observed previously 
in the simulation of systems with COP and without topological defects[7]. 



16 



Finally, in agreement with Rojas and Bray[12] we find that the peak of F[x) 
is well fitted by the quartic exponential form appearing in the BH analytical 
solution. 



5 - Conclusions 



The main motivation for this paper was to investigate in detail the onset 
of standard scahng in the BH model for phase ordering kinetics with COP 
and finite N. We have done this by a comparative study of the numerical 
solution of the model with NCOP and with COP. In both cases eventually 
there is standard scaling, but the difference is much more profound than just 
the value of the growth exponent (z = 2 for NCOP and z = 4 for COP) when 
the whole development of the dynamics is taken into account. As probes for 
scaling we have used Y{t) and (p{x). Parameters of the quench are the initial 
condition A and the number of components N of the order parameter. 

The picture for NCOP is the following. Starting from a uniform initial 
condition C{k,t — 0) = A, after a short transient of duration to during which 
information on the initial condition is lost, the dynamics of standard scaling 
sets in, with ordered regions growing like L{t) ~ t^/^, (p{x) = d and the 
scahng form (2.6) obeyed. The only place where there remains a detectable 
transient dependence on the initial condition A for longer times than to is in 
the behaviour of Y{t) which displays a very slow approach to the constant 
asymptotic behaviour. This means that in the scaling ansatz for R(t) there 
is a slow correction with a small amplitude. This pattern of behaviour for 
NCOP is the same for any value of N, including N — oo. 
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By contrast, with COP the way the system eventually reaches standard 
scaling is more complicated and depends on N due to the existence of the 
two characteristic times to and t*. Only if t* ~ to there is no observable 
difference between COP and NCOP. Instead, if t* is sufficiently larger than 
to the system displays multiscaling in between to and t*, before the standard 
scaling regime is reached. In this sense multiscaling is not only a feature 
of the special case N — oo, but is relevant also for systems with finite N. 
The existence of a connection between multiscaling and standard scahng is 
expected after recognizing that these are the two asymptotic features of a 
crossover process. More specifically, let us consider a general scaling form 
containing both regimes 

C{k,t,N)^L'^(^^J^{x,^) (5.1) 

with X — k/km, a an index to be determined and T a function with the 
limiting behaviors 

-rr il^-/ 1 if7V/L"»l 

~ \ A{^)'^(^) if A^/L" < 1 ^^-^^ 

where A and a{x) also must be determined. The above form clearly yields 
multiscaling if the limit — >^ oo is taken. If instead N is kept finite and L{t) 
gets large one has 

C{k, t, N) = ^L'^V'{a;)-aaW7V"W (5.3) 

and imposing dip{x) — aa{x) — d one finds standard scaling 

C{k,t,N) = AL'^e-^(^-'^(^»i"^ = AL'^e-^(^'-i)'i^^ (5.4) 



exactly with the BH scaling function revealing the deep connection between 
multiscaling and standard scaling as the multiscaling spectrum ip{x) dictates 
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the form of the scahng function in the standard scahng regime. This mul- 
tiscahng to standard scahng crossover in principle could be observed also 
in systems with realistic values of N by making t* large enough exploiting 
the dependence of t* on A. In this respect it might be interesting to check 
this hypothesis on the simulations of ref.s [5,6,7] performed with values of A 
making t* sufficiently large. 

Finally, the finding of exponential tails in the scaling functions is quite 
interesting and CDS simulations are under way in order to check on the 
existence of these tails in systems with N > d + 1. 
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Figure Captions 

Fig.l - lnY(t) for NCOP with N = oo and different values of A displaying 
tfie approach to ln(l/4) = -1.38629 for d = 1 and to ln(l/2) = -0.69314 for 
d — 2 as predicted by (2.15). 

Fig.2 - \nY{t) for COP with = oo and different values of A. 

Fig. 3 - Time evolution of (p{x) for COP with N — oo, d — 1 and different 

values of A. 

Fig.4 - Late stage (200000 <t < 500000) multiscaling behaviour of (p{x) for 
COP with N = oo revealing independence from the initial condition. 
Fig.5 - InF(t) for NCOP with = 10, 10^ 10^ The A dependent transient 
preceding the asymptotic behaviour is negligible in this scale for d = 1. 
Fig. 6 - Plot of the scaling function for NCOP demonstrating exponential 
decay in the tails. 

Fig. 7 - laY{t) for COP with N — 300 displaying, in the preasymptotic 
regime, the switch from standard scaling to multiscaling with increasing A. 
Fig. 8 - laY(t) for COP with A = 0.01 displaying, in the preasymptotic 
regime, the switch from standard scaling to multiscaling with increasing N. 
Fig. 9 - The evolution of Lp{x) for COP with A = 0.01, li = 2 and different 
values of A^. 

Fig. 10 - Plot of the scaling function for COP demonstrating exponential 
decay in the tails. Computations have been carried out with A = 0.01. 
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